clear; clc; close all;
% 在 c=0.2, x0=0.1, q0=0.1 条件下的吸引子图，
% HCA (Hidden Chaotic Attractor): 指系统出现的隐藏混沌吸引子
% HHCA (Hidden HyperChaotic Attractor): 指系统出现的隐藏超混沌吸引子。

% 模型参数
a = 0.1;
b = 0.1;
c = 0.2;               % 固定 c=0.2
k_values = [1.75, 1.85];% k=1.75 和 k=1.85
x0 = 0.1;              % 初始 x
q0 = 0.1;              % 初始 q
transient = 100;      % 舍弃暂态步数
iterations = 20000;    % 记录阶段迭代步数

% 预分配存储
x_rec = zeros(iterations, length(k_values));
q_rec = zeros(iterations, length(k_values));

% 迭代并记录
for idx = 1:length(k_values)
    k = k_values(idx);
    % 初始化当前状态
    x_curr = x0;
    q_curr = q0;
    
    % 暂态阶段
    for n = 1:transient
        [x_curr, q_curr] = mptm(x_curr, q_curr, a, b, c, k);
    end
    
    % 正式记录阶段
    for n = 1:iterations
        [x_curr, q_curr] = mptm(x_curr, q_curr, a, b, c, k);
        x_rec(n, idx) = x_curr;
        q_rec(n, idx) = q_curr;
    end
    
    % 绘图
    subplot(1, 2, idx);
    plot(x_rec(:, idx), q_rec(:, idx), '.', 'MarkerSize', 1, 'Color', 'g');
    xlabel('x');
    ylabel('q');
    if idx == 1
        title('(a) k=1.75, HCA');
    else
        title('(b) k=1.85, HHCA');
    end
    axis tight;
    grid on;
end
